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S-h ' Abstract. We consider the problem of electron energy states related to strongly 

^ ■ localized potential of a single impurity in graphene. Our model simulates the effect of 

impurity atom substituting the atom of carbon, on the energy spectrum of electrons 
r^ near the Dirac point. We take into account the internal spin-orbit interaction, which 

1 , can modify the structure of electron bands at very small neighborhood of the Dirac 

point, leading to the energy gap. This makes possible the occurrence of additional 
f^ impurity states in the vicinity of the gap. 
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1. Introduction 

Graphene attracted a lot of attention recently due to very unusual properties of electron 
energy spectrum and transport properties, including both the transport of electrons and 
phonons [U El EJ EJ IS]. The most striking properties of graphene are related to the 
energy spectrum near the Dirac points, where this spectrum is linear as a function of 
momentum, and the Hamiltonian of free electrons can be described by the relativistic 
two-dimensional Dirac model [6]. 

Naturally, the impurities and defects can strongly affect the energy spectrum of 
graphene. Especially important is the effect of impurities on the spectrum near the 
Dirac point. The impurity states and the corresponding variation of the electron density 
of states have been already discussed in several papers [3 El El [10] without taking 
into account the spin-orbit (SO) interaction. It was found that the localized impurity 
potential gives the resonant states in the spectrum of graphene. They can be located 
near the Dirac point in the case of relatively strong impurity potential, and this is 
quite unusual for the semiconductor physics. In the case of carbon vacancy, there 
appears a local energy level at the Dirac point, with E = 0. In all of these works, the 
main attention has been paid to a finite density of impurities and defects leading to 
modification of the density of states in graphene. 

In this paper we mostly concentrate on the problem of single impurity taking into 
account the internal SO interaction. The SO interaction opens a gap in the electron 
energy spectrum [IT]. However, it was found that the magnitude of SO-induced gap is 
very small in graphene [121 HSl HH US] , an d therefore it would be very difficult to observe 
this gap experimentally. Nevertheless, the problem exists: how this small gap would 
affect the behavior of impurity states in the vicinity of Dirac point? 

We demonstrate that the SO gap induces appearance of additional impurity states 
corresponding to very weak impurity potential. As a result, the SO gap in graphene can 
be experimentally unobservable due to a large number of adsorbed light atoms at the 
graphene surface creating the impurity states in the gap. 

2. Model 

We use the following Hamiltonian, which describes electrons near the Dirac point K, 
with the intrinsic SO interaction |TTJ 




h* = r ":- , (i) 



where A is the band splitting related to the SO interaction, v is the velocity parameter, 
and we denote k± = k x ±ik y . The matrix form of (0Q) is due to the choice of wavefunction 
basis corresponding to different sublattices in the lattice of graphene. The basis functions 
of Hamiltonian dTJ) are 

|k(l,2)a)= £ e «< k » +k >"V.(r) k) , (2) 

se(A,B) 
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where ip s { r ) is the tight-binding electron state at site s belonging to sublattice A or B, 
k is the wave vector corresponding to the chosen Dirac point, and a =f, \, refers to spin 
up and down states, respectively. The eigenvalues of Hamiltonian (CQ) are Eki,2 = ±£fc, 
where e^ = (A 2 + v 2 k 2 ) , so that the value of 2| A| is the energy gap. The Hamiltonian 
for the other Dirac point K' differs from (CQ) by the opposite sign in the diagonal terms. 
Thus, the results for the point K' are can be found by using calculations with the 
Hamiltonian ([1]) and reverting the sign of A. 

Since the SO interaction in Equation (TjQ) does not mix the spins, one can consider 
separately spin up and down channels. For the spin up electrons the Hamiltonian can 
be presented as it gives the Hamiltonian for up-spin electrons 

# kt = T Z A + VT-k, (3) 

where r, are the Pauli matrices acting in the space of sublattices A and B. For the 
down-spin Hamiltonian, the sign of A in Equation ([3]) is opposite. We consider first the 
spin up Hamiltonian (J3D- 

3. Nonmagnetic impurity 

Let us consider the impurity state in the case of a single impurity, described by 
the perturbation localized in one of the sublattices. In the continuous model under 
consideration it corresponds to the perturbation at r = in sublattice A 

Mr) = ( V °T ] ) ' <«> 

The matrix of perturbation (J3J) in the basis functions of Hamiltonian H^ is 

fc*t =( V n n ) = *W- (5) 



The effect of perturbation Vkk' t o n the energy spectrum in all orders of magnitude can 
be described using the T-matrix method [IB] . In the general case, the equation for the 
T-matrix is 

TMe) = VW + £tW6 kl (e)t klk ,(e), (6) 

ki 

where G k (e) — [e — H^) is the Green's function of Hamiltonian H^. Using ([3]) we 



find the Green's function for spin up electrons 

G ^( e ) = 7^—2 • ( 7 ) 

In the following we assume that the energy parameter includes a small imaginary part, 
e — > e + iSsgne, which corresponds to the choice of retarded Green's function. Using 
(ED and (ED we find 

-i 



T t (e) 



i-yot£&t( £ ) 



v ( 



ot 



V ot F t {e] 



-i 



Vot, (8) 
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and we have to calculate 

k 

e + i5sgu£ + r z A , v 2 k 2 , + A 2 - e 2 - 2i\e\8 
— In 



Attv 2 



A 2 



2i\e\5 



(9) 



where the upper limit (cutoff) of integration over momentum k m is introduced. It 
corresponds to the region of linearity of the spectrum near the Dirac point. Assuming 

v 2 k 2 > I A 2 — e 2 \ we obtain 



F,(e) 



e + i^sgne + t z A 
Attv 2 



2h,2 



In 



v 2 k. 



(A 2 - e 2 ) 2 + Ae 2 8 2 



1/2 



+ i(<pi -<f 2) 



(10) 



where fi^ are the angles related to the phase of complex function F(e), which can be 
made analytical in the whole complex plane of e after proper choice of cuts in this plane. 
We assume that te cut is made along the real axis from — | A | to + 1 A | . Then the phases 
can be chosen 



<Pi 



f2 



— tan 



-1 2|e|* 
v 2 k 2 



7r + tan x 



2\e\S 



A 2 



tan" 



2\e\S 



\e 2 - A 5 



£ 2 >A 2 



e 2 < A 2 . 



:u) 



(12) 



The real part of Equation 



ReF t (e) 



e + r z A 
Attv 2 



In 



27.2 

■in 



v 2 k 



1/2 • 



(A 2 - e 2 y + Ae 2 5 2 
In the limit of S — > (corresponding to small impurity density) 

ReF t (e) 



e + r z A v 2 k 2 m 
~ In — 



ImF t (e) 



Attv 2 

e + r z A 
Av 2 '' 



A 2 -e 2 \ 



e 2 >A 2 



(13) 



(14) 



(15) 



0, e 2 < A 2 , 

The functions Rei r n j -j-(e) and ImFn^(e) are presented in figure [Hand figured 

The matrices Vot anc i -^t(e) are both diagonal. Therefore, T{e) calculated from 
Equation ([8]) is diagonal, too 

Vot 



?M 



diag 



,0 



v 1 - Vot Re ^ii,t( e ) ~ iV ot Im ^ii,t( e ) 
The location of impurity level is determined by the pole of T-matrix 



1-VWReiWe 



0. 



(16) 



(17) 
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Figure 1. The dependence ReFn(e) for spin-up and spin down states, related to the 
Dirac point K in the case of non-magnetic impurity. 



5 


X 


0" 


















* 








i 




: 


3 








•=r> 




^~^"^' 


*> 1 

(H 




....--" 










" 


-3 







-0.01 



-0.005 



0.005 



0.01 



£c 



Figure 2. The imaginary part ImFn(e) for spin-up and spin down states, related to 
the point K in the case of non-magnetic impurity. 



The dependence ReF ll:1 -(e) presented in figured] can be used for the graphical solution 
of Equation ( tTTjl . Using ( TT4l) we find 



l + ^l+Aln *%. 



'Of 



4irv 2 |A 2 -e 2 | 
Thus, the equation for the impurity level is 

A.'KV 2 

e t = -A 



0. 



v 2 k 2 



(19) 



A 2 -4 



This equation has several solutions for the same potential Vbf . 
The Hamiltonian describing the spin down states 



#i 



14 



-7\A + wt ■ k, 



(20) 
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differs from ([2]) only by the sign of gap parameter A. In the case of nonmagnetic 
impurity the perturbation Vj,(r) = Vj-(r). Performing the same calculations as before 
with the substitution A — > —A we find the equation for impurity level corresponding 
to the spin-down state 

Anv 2 



A 



21.2 



Vnln 



v 2 k 



(21) 



'i 



The spin up and down states corresponding to solution of Equation ([19]) and Equation 
( T2~Tj) . respectively, are split in energy. The magnitude of splitting is 2|A|, which of 
course is very small as was discussed in the Introduction. Nevertheless, considering the 
impurity states in the model with one Dirac point we come to a weakly magnetized 
state at the nonmagnetic impurity. This nonequivalence of spin up and down states is 
exactly compensated by the states related to another Dirac point, K', for which the 
Hamiltonians of spin up and down states differ by the sign of A from those in Equation 
([3]) and Equation (|2"T1) [Tlj. Thus, the magnetization of the localized state is absent if 
we take into account both nonequivalent Dirac points. 
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Figure 3. The location of impurity levels as a function the impurity strength 
parameter Vb (a) and the schematic presentation of the impurity levels (b) in the 
case of non-magnetic impurity. 



The numerical solutions of Equation (fT9"|) and Equation 02 1|) describing the states 
related to the Dirac point K, as well as the corresponding solutions related to the 
Dirac point K' are presented in Figure EJa. The schematic representation of the 
levels is presented in Figure [3],b. All the levels are spin degenerate due to the 
overlapping solutions of equations, related to the nonequivalent Dirac points. We used 
the parameters with much larger SO gap A to visualize better the character of solutions. 
If the state is located within the gap, the impurity level is discrete. For the level with 
energy \e^\ > A, it is a resonant state of width lm.F^(e^). It should be noted that 
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ImF^(e) is not the density of states in graphene because Equation Q does not include 
the trace over sublattices. 

It should be noted that both impurity states with energies e-j- and £4 are mostly 
localized at the site of A sublattice in accordance with the assumption that the impurity 
potential fl3]) is located on the A-site. However, the real spread of the wavefunction can 
be much larger than the distance between nearest A and B sites. 

4. Magnetic impurity 

In the case of magnetic impurity we choose perturbation, which is different in sign for 
the spin up and down electrons, VWj, = ~~ ^kk't> an d ^kk't is described by Equation ((5]). 
The resulting impurity levels related to the Dirac points K and K' are to be found from 
the following equations 

W . (22) 



: t4 



TAT 



21.2 



Vnln 



v 2 k. 






±A T 




2ju2 



Vnln 



v k 



(23) 



A 2 



't4 



There is the spin splitting of the states related to the K point, which does not vanish 
as A — y 0, and the resulting magnetization of impurity states is not compensated by 
another point K' . Corresponding numerical solutions of equations ( 122|) and (|23|) are 
presented in Figure HI 
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Figure 4. The same dependences as in Figure [3] in the case of magnetic impurity. 
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5. Nonmagnetic impurity with strong SO interaction 

Using the same formalism, one can consider the impurity, which locally enhances the 
internal SO interaction. In this case the corresponding perturbation for K point has 
the form 

where A is the local After substitution of Vq — > A , the solution for the K point does 
not differ from the case of magnetic impurity considered above. However, for the point 
K' we have to change the sign of Ao- As a result we obtain 

Airv 2 



't4 



J? 



tat 



An In 



±A± 




2 7.2 



An In 



v k, 



(25) 



(26) 



A 2 



: t,4- 



In this case the spin splitting of the states related to the point K does not vanish as 
A — > but the resulting local magnetization is absent due to the states from point K' . 
The solutions of equations ( 1251) and ( 1261) are presented in figure 
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Figure 5. The same dependences as in Figure |3] in the case of spin-orbit impurity. 



6. Wave function of the localized impurity state 

Using the Schrodinger equation for the wave function, one can find the equation for 
impurity state [IS] 

^i(r) = fd 2 r'G(r,r';e)V(r')Mr'), (27) 
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where the Green function of free electrons in coordinate representation obeys 

(e-H r )G(r,r'-£) = S(r-r'). 
For the impurity located at r = we obtain from ( 1271) 

^(r) = Gf(r,0;e)f(e)^ o (0), 
where ipo(r) is the eigenfunction of unperturbed Hamiltonian (OQ) 
(e-H r ) Vo(r) = 0. 
The Green function for Hamiltonian ([3]) has been calculated before [TTJ 

Ar z 



(28) 
(29) 
(30) 



G t (r,r';e) 



1(8 



Av 2 
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A 2 



+ r 



Ve 
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Av 2 \v — v 



7^ 



(i) 



|r-rV^ 2 - A 2 



v 



Thus, for the spin-up impurity state with energy e-j- 

2 



i^twr 



a, 



ID 





"t 



m 



(i) 



A s 



-2r/R 



(31) 



(32) 



«/ 



where HJP(z) are the Hankel functions [18] and i?ot 
characteristic radius of the impurity wavefunction. 

The function ^o( r ) in Equation ( 1291) can be calculated using ( 130 
bispinor components ipQ (r) = ((^ T (r), x T ( r ) ) we find the relation 



A 5 



1/2 



is the 



Denoting the 



X(r) 



—iv(d x + i<9 y 



</?(r) 



A + e 
The equation for the function tp^ with the Hamiltonian 

reads 

/ d 2 Id Id 2 



dr 2 r dr r 2 da 2 



«t ^t( r > «) = ° 5 



which has the following solutions decaying at large r 

p t (r,a) = J M« t r)e ima , 
where «^ = (A 2 — si) /v 2 , K m (z) is the modified Bessel function and m E Z. Then the 



(33) 
in polar coordinates (r, a) 

(34) 
(35) 



(36) 



(37) 



corresponding ^-component can be found from ( 133]) and ([35 
X,(r,a) = -^K m+1 (K t r)e^«. 

l\ + £-)- 

Thus, up to the normalization we have 

It should be noted that for ^o(O) in Equation (]2"9"]) we have to use the cutoff at small 
distance r c ~ k^ ax ^> a /V, where a is the lattice constant and k max is the upper limit 
for the Dirac model in graphene. It means that using the Dirac model we can find the 
impurity wave function at distances much larger than Oo- On the other hand, as we see 
from (]32l . the radius of the impurity states Rq near the Dirac point is much larger than 
this limit. 
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7. Conclusions 

We have calculated the energies and wave functions of impurity states near the Dirac 
points in graphene taking into account the SO interaction. The calculations show 
the SO-induced spin splitting of these states. The existence of two nonequivalent 
Dirac points in the Brillouin zone leads to the spin degeneracy of the states with 
different spin. It should be noted that in principle the valley degeneracy can be 
broken by inhomogeneous deformations, which would result in the appearance of local 
magnetization. 
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